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For a weakly coupled quantum field at high temperature the classical approximation offers a possibility to 
gain insight into nonperturbative real-time dynamics. I use this to present a nonperturbative approach to the 
computation of spectral functions in real time. Results are shown for a scalar field in 2 + 1 dimensions. 



1. In thermal equilibrium spectral functions 
contain all available information since other real- 
(and imaginary-) time correlators are related to 
them via the KMS condition. A nonperturbative 
computation of spectral functions seems therefore 
desirable; nonperturbative in order to: 

i) deal with a truly nonperturbative the- 
ory, such as the infrared sector of high- 
temperature nonabelian gauge fields, 

ii) investigate the utility of (resummed) per- 
turbative ideas, since weak coupling + high 
temperature ^ simple. 

For the calculation of static quantities, such 
as free energies, phase diagrams and screening 
masses, the euclidean lattice formulation of finite- 
temperature quantum fields offers a first-principle 
approach. For intrinsically real-time correlators, 
such as the spectral function, no first-principle 
nonperturbative formulation exists. In the last 
O(10) years the classical approximation at high 
temperature and weak coupling has provided 
much insight into real-time dynamics, see |j],|2| for 
recent reviews. In this talk I discuss how spectral 
functions can be computed in the classical ap- 
proximation and what can be learned from them. 
Details can be found in (D • 

2. The spectral function for a bosonic operator 
O can be defined as the thermal expectation value 



p{x-y)=i{[0{x) > 0\y)]. 



(1) 



Other two-point functions can be expressed in 
terms of it using the KMS condition, leading 
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to the appearance of the Bose distribution func- 
tion n{uj) = l/[e"/ T — 1]. For example, the eu- 
clidean correlator D(t, x) = (0(t, x)0' f (0, 0)}e, 
with imaginary time r € [0, 1/T], is related to 
the spectral function via the integral equation 



with the kernel 
K(t, to) = tl(uj) 



-K(r,uj)p(u;,p), 



-u(t-1/T) 



(2) 



(3) 



obeying K(t,l>) = -K{t,-w) = K(l/T - r,w). 
The extraction of the spectral function from the 
euclidean lattice correlator is a highly nontriv- 
ial inversion problem, currently tackled using 
the Maximal Entropy Method 0. In real time 
the spectral function would be accessible directly 
without the need to solve an integral equation. 

The classical approximation provides a nonper- 
turbative computational scheme in real time. In 
the classical limit —i times the commutator is re- 
placed with the Poisson bracket and the classical 
spectral function reads 



Pci(x-y) = -({0(x),0^y)}) ch 



(4) 



The time evolution is determined by the classical 
equations of motion, evolving from some initial 
condition. The brackets (-) c i denote a Boltzmann- 
weighted average over these initial conditions. 

In perturbation theory the relation between 
classical and quantum physics at finite temper- 
ature can be seen clearly. Consider a self- 
interacting scalar field with coupling constant A. 
The effective expansion parameter is h\n(huj) 
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Figure 1. Spectral functions p c \{oJ, 0) for various 
temperatures T. Fits to a Breit-Wigner function 
are shown with dotted lines. 
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Figure 2. As in Fig. 1 at T/m = 7.2. The actual 
data points are presented with squares. 



(with the h dependence indicated explicitly). 
As a result, in the classical limit each loop 
contributes with a proportionality factor XT/co. 
Note that this implies that the infrared sector 
(Ticu/T <C 1) is essentially classical but that 
in the ultraviolet problems are encountered (the 
Rayleigh- Jeans divergence, see below). 

For actual calculations, the Poisson bracket 
in Eq. (4) appears difficult to use numerically, 
but in equilibrium the KMS condition comes to 
the rescue. In the classical limit the KMS rela- 
tion relates the statistical correlator S(x — y) = 
(0(x)0^ (y)) c i and the spectral function in mo- 
mentum space as 

T 

- i-p c i(u),p) = S(u),p). (5) 

LO 

or equivalently in real-space as 

p d (t,x) = ~d t S(t,x), (6) 

The latter relation, valid for arbitrary O, allows 
an easy calculation of classical spectral functions. 

3. To illustrate the method we consider a simple 
real scalar field with a A^ 4 /4!-interaction in the 
symmetric phase in 2 + 1 dimensions (note that in 
a classical theory the dimcnsionless coupling A/to 
can be taken equal to 1 without loss of generality) . 
We focus on the one-particle spectral function. 



According to the classical KMS condition we have 
to compute 

p cl (t,x) = -<{0(t,x),0(O,O)}) c i 

= -i<7r(i,x)<K0,0)) cl , (7) 

where 7r = dt<fi is the canonical momentum. The 
right-hand-side is a simple correlation function 
that can be computed numerically without any 
problems ||. The theory is defined on a spatial 
lattice of N x N sites and lattice spacing a (we 
used N = 128, ma = 0.2, and periodic boundary 
conditions). The classical equations of motion are 
solved with a leapfrog algorithm with time step 
ao/a = 0.1. Thermal initial configurations are 
generated with the Kramers equation algorithm 
(the plots shown are obtained with 2000 indepen- 
dently thermalized initial configurations). On the 
lattice the spectral function (7) is symmetrized: 

Pci,iat(£,x) = 

~\ (*(t + \a , x)i [<X0, 0) + 0(a o , 0)Yj . (8) 

The simulations give the spectral function as a 
correlation function in real time. A sine trans- 
form (using symmetry properties of the spectral 
function under time reflection and complex con- 
jugation) yields the desired result as a function of 
frequency. In Fig. 1 the spectral function at zero 
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Figure 3. Temperature dependence of the mass 
M (circles) and width T (squares, multiplied by 
25 for clarity), obtained from fits to a Breit- 
Wigner function. Errors are determined with a 
jackknifc analysis. The resummed perturbative 
predictions are shown with lines. 



momentum is shown for various temperatures and 
in Fig. 2 a magnification with the actual data 
points indicated is presented. The maximal (real) 
time i max used in the analysis constrains the res- 
olution in frequency space to Au> = ir/t max , but 
this poses no problem (Aw/m ~ 0.01 in Fig. 2). 

4. It is now clear that spectral functions can be 
computed nonperturbatively in classical thermal 
field theory. The question is what we may learn 
for hot quantum fields. 

As alluded to above, the setup of perturbation 
theory in the quantum and classical case is sim- 
ilar and the results may be used to address the 
applicability of (resummed) perturbation theory. 
In the weak-coupling regime the spectral func- 
tion is dominated by the plasmon. A simple way 
to parametrize it (ignoring the proper analytical 
structure) is by a Breit-Wigner spectral function. 
At zero momentum it reads 

PBW{UJ > ~ (^-M^+^' (9) 

Fits of the data to a Breit-Wigner function 
are shown in Figs. 1,2, and yield nonperturba- 
tively determined values of the plasmon mass and 
width. These can be compared with the pertur- 



bative predictions. At lowest order the mass pa- 
rameter M is determined from the classical limit 
of the standard one-loop gap equation. The lead- 
ing contribution to the width T comes from the 
two-loop setting-sun diagram and on dimensional 
grounds it is T = c\ 2 T 2 /M 3 , with a small coeffi- 
cient c = (3 - 2^/(32^) §. A comparison be- 
tween the perturbatively and nonperturbatively 
determined values of the effective mass and width 
is shown in Fig. 3. A nice agreement can be seen, 
indicating that for this range of parameters per- 
turbation theory is reliable. 

5. It is straightforward to extend the calcula- 
tion to more complicated spectral functions in 
scalar and gauge theories (for lattice fermions 
in real time, see ||). Transport coefficients 
can be defined from the zero-momentum and 
zero-frequency limit of equilibrium spectral func- 
tions of appropriate composite operators. Simple 
power counting shows that these quantities are 
dominated by hard (~ T) momenta. In a clas- 
sical calculation they will therefore be sensitive 
to the lattice regulator. In fact, such dependence 
has already been encountered since the effective 
mass parameter M depends (logarithmically in 
2 + 1 dimensions) on the lattice cutoff. However, 
as can be seen from Fig. 3, this does not automat- 
ically imply that the classical findings cannot be 
used: analytical perturbative and numerical non- 
perturbative calculations can still be compared, 
provided the role of the lattice regulator is incor- 
porated properly. 
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